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Abstract: We demonstrate how to realize within supergravity a novel chaotic-type infla¬ 
tionary scenario driven by the radial parts of a conjugate pair of Higgs superfields causing 
the spontaneous breaking of a grand unified gauge symmetry at a scale assuming the value 
of the supersymmetric grand unification scale. The super potential is uniquely determined at 
the renormalizable level by the gauge symmetry and a continuous R symmetry. We select 
two types of Kahler potentials, which respect these symmetries as well as an approximate 
shift symmetry. In particular, they include in a logarithm a dominant shift-symmetric term 
proportional to a parameter c_ together with a small term violating this symmetry and 
characterized by a parameter c+. In both cases, imposing a lower bound on c_, inflation 
can be attained with subplanckian values of the original inflaton, while the corresponding 
effective theory respects perturbative unitarity for = c+/c_ < 1. These inflationary models 
do not lead to overproduction of cosmic defects, are largely independent of the one-loop 
radiative corrections and accommodate, for natural values of r±, observable gravitational 
waves consistently with all the current observational data. The inflaton mass is mostly 
confined in the range (3.7 — 8.1) x 10^° GeV. 
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1. Introduction 

The announcement of the recent joint analysis of the Bicep2/&c^ Array and Planck data [1,2] 
confirms earlier attempts [3,4] which stressed the impact of the dust foreground on the observations on 
the B-mode in the polarization of the cosmic microwave background radiation at large angular scales. 
As a consequence, the predicted tensor-to-scalar ratio r in inflationary models must be smaller than the 
one initially claimed in Ref. [5]. However, the present data not only leave open the possibility for a 
sizable value of r, but also seem to favor values of r of order 0.01. Indeed, it has been reported that 

r = 0.048l°;0| (1.1) 

at 68% confidence level (c.l.). This fact motivates us to explore the question whether realistic su¬ 
persymmetric (SUSY) inflation models can accommodate such values of r - for similar attempts see 
Refs. [6-10]. 

One elegant SUSY model which can nicely combine inflation with the Higgs mechanism of the 
symmetry breaking is the model of Higgs inflation (HI). This is an inflationary model of the chaotic 
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type, where a Higgs field plays the role of inflaton before its trapping in the vacuum. It has been 
shown that HI in the framework of supergravity (SUGRA) can be implemented by imposing [11-15] a 
convenient shift symmetry on the Kahler potential or invoking [16-18] a logarithmic Kahler potential 
with a real subdominant kinetic part and a dominant holomorphic (and anti-holomorphic) part, which 
plays the role of a non-minimal coupling to the Ricci scalar curvature [19, 20]. Inspired by these 
efforts, we present here a ‘hybrid’ scenario, where a conjugate pair of Higgs superfields is involved in 
the logarithmic part of the Kahler potential K, which respects a shift symmetry with a tiny violation 
- cf. Refs. [10, 21, 22]. The shift-symmetry-preserving part of K influences fhe amplitude of the 
canonically normalized inflaton, which becomes much larger than the original inflaton field appearing 
in the superpotential and the Kahler potential - cf. Ref. [8]. Therefore, HI can be implemented even 
with subplanckian values of the original inflaton field, keeping corrections from higher order terms 
harmless. On the other hand, the resulting inflationary potential does not depend directly on the strength 
of the shift-symmetry-preserving part of K and, thus, is not flattened drastically as in the original 
scenario of non-minimal HI [16-18], but just adequately by the shift-symmetry-violating part of K, as 
in the recently proposed models of kinetically modified non-minimal inflation [10]. Moreover, invoking 
deviations from the prefactors —3 or —2 of the logarithms in the proposed Kahler potentials, we succeed 
to enhance the resulting values of r - cf. Refs. [6, 8,9]. We also analyze the impact of the one-loop 
radiative corrections (RCs) [23] on our results and find that these can be kept under control provided 
that the relevant renormalization scale is conveniently chosen [24]. We, finally, show that the ultraviolet 
(UV) cut-off scale [25-27] in these models coincides with the Planck scale and so concerns regarding 
their naturalness can be safely eluded. 

We exemplify our proposal in the context of a grand unified theory (GUT) model based on the 
gauge group Gb-l = Gsm x U{1)b-l - where Gsm = 5 '( 7 ( 3 )c x SU{2)i^ x U{1)y is the gauge 
group of the standard model and B and L denote the baryon and lepton number, respectively. Actually, 
this is a minimal extension of the minimal supersymmetric standard model (MSSM) which is obtained 
by promoting the already existing U{1)b-l global symmetry to a local one. As a consequence, the 
presence of right-handed neutrinos vf is necessary in order to cancel the B — L gauge anomaly. The 
Higgs fields which cause fhe spontaneous breaking of the Gb-l symmetry to Gsm can naturally play 
the role of inflaton. This breaking provides large Majorana masses to the uf’s, which then generate the 
tiny neutrino masses via the seesaw mechanism. Furthermore, the out-of-equilibrium decay of the uf’s 
provides us with a robust baryogenesis scenario via non-thermal leptogenesis [28]. 

It is worth emphasizing that U{1) b-l is already spontaneously broken during HI through the non¬ 
zero values acquired by the relevant Higgs fields. Consequently, HI is not followed by the production 
of cosmic strings and, therefore, no extra restrictions [29] on the model parameters have to be imposed, 
in contrast to the case of the standard F-term hybrid inflation (FHI) [30] models, which share the same 
superpotential with our models. In the standard FHI models, the GUT gauge symmetry is unbroken 
during inflation since the Higgs superfields are confined to zero and the inflaton is a gauge singlet. 
The spontaneous breaking of the GUT gauge symmetry takes place at the end of FHI, where the Higgs 
fields acquire non-zero values. Topological defects are, thus, copiously formed if they are predicted by 
the symmetry breaking. In our present scheme, this same gauge singlet superfield is sfabilized at zero 
during and after HI. We consider two possible embeddings of the gauge singlet superfield in K with its 
kinetic terms included or not included in the logarithm together with those of the inflaton. 

The superpotential and Kahler potentials of our models are presented in Sec. 2. In Sec. 3, we de- 
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scribe the inflationary potential at tree level and after including the one-loop RCs, whereas, in Sec. 4, we 
derive the inflationary observables and confront them with observations. We then provide an analysis 
of the UV behavior of these models in Sec. 5. Our conclusions are summarized in Sec. 6. Throughout 
the paper, we use units where the reduced Planck scale mp = 2.433 x 10^® GeV is set equal to unity, 
subscripts of the type , x denote derivation with respect to (w.r.t.) the field x (e.g. = d^Fjdx^), 

and charge conjugation is denoted by a star. 

2. Modeling Shift Symmetry for Higgs Inflation 

We will now explain how a conveniently modified shift symmetry can be used in order to imple¬ 
ment HI based on the F-term SUGRA potential. The structure of the superpotential is presented in 
Sec. 2.1, whereas the relevant Kahler potential is given in Sec. 2.2. Finally, in Sec. 2.3, we derive the 
corresponding frame function. 

2.1 The superpotential 

We focus on a minimal extension of the MSSM based on 
the gauge group Gb-l, which can be broken down to Gsm at a 
scale close to the SUSY GUT scale Mqut through the vacuum 
expectation values acquired by a conjugate pair of left-handed 
Higgs superfields and l> charged oppositely under U{1)b-l Table 1: Charge assignments of the 
- see Table 1. The part of the superpotential W which is relevant superhelds. 
for inflation is given by [30] 

w = xs , ( 2 . 1 ) 

where 5 is a gauge singlet superfield, A a dimensionless parameter, and M a mass scale of order Mgut- 
This superpotential is the most general renormalizable superpotential which respects an R symmetry 
U{1)r - see Table 1 - in addition to the aforementioned Gb-l- The R symmetry guarantees the 
linearity of the superpotential w.r.t. the gauge singlet superfield S. This fact is helpful both for the 
realization of HI and the determination of the SUSY vacuum. 

To verify that W leads to the breaking of Gb-l down to Gsm> we minimize the SUSY limit 
Tsusy of the SUGRA scalar potential derived from the superpotential in Eq. (2.1) and the common 
SUSY limit of the Kahler potentials in Eqs. (2.8) and (2.9) - see below. The potential VsusY? which 
includes contributions from E- and D-terms, turns out to be 

^SUSY = 1.$ - ^ (|$|2 + |ch|2) + - Nr±f (|l>|2 - , 

(2.2a) 

where the complex scalar components of the various superfields are denoted by the same superfield 
symbol, g is the unified gauge coupling constant, and the remaining parameters {N, c_, r±) are defined 
in Secs. 2.2 and 2.3. Erom the last equation, we find that the SUSY vacuum lies along the D-flat 
direction |$| = |<I>| with 

(S) = 0 and |($)| = |(|.)| = M/2. (2.2b) 

Although (<I>) and ($) break spontaneously U{1)b-l, no cosmic strings are produced at the end of 
inflation, since this symmetry is already broken during HI. Needless to say that contributions from 
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the soft SUSY breaking terms ean be safely negleeted sinee the eorresponding mass seale is mueh 
smaller than M. Let us emphasize, however, that soft SUSY breaking effeets break U{l)ji explieitly 
to a diserete subgroup. Usually, the eombination of the latter with the fermion parity yields [31] 
the well-known i?-parity of MSSM, whieh guarantees the stability of the lightest SUSY partiele and, 
therefore, provides a well-motivated eold dark matter eandidate. 

2.2 The Kahler Potential 

The superpotential W in Eq. (2.1) eould give rise to HI driven by the real field (/> defined in fhe 
sfandard paramefrizafion 

^ =-^e*®cos0<i>, ^ = -^e*^sin0# wifh 0 < 0$ < ^ and S = ^ (2.3) 

\/2 \/2 2 y 2 

provided fhaf we eonfine ourselves fo fhe field eonfigurafion 

s = s = 9 = 9 = 0 and 0$ = vr/d. (2.4) 

Nofe fhaf fhe lasf equalify ensures fhe D-flafness of fhe pofenfial. Indeed, along fhis frajeefory, VsusY 
in Eq. (2.2a) reduees fo fhe well-known E-ferm pofenfial whieh is quarfie w.r.f. (p. This eonsfrue- 
fion, fhough, ean beeome meaningful only if if ean be sueeessfully embedded in SUGRA, due fo fhe 
fransplanekian values of p whieh may be needed - see below. 

To fhis end and following similar works [11,13,15], we require fhaf fhe Kahler pofenfial is eonsis- 
fenf wifh fhe shiff symmefry 


$ ^ ^ ^ ^ + C*, S (2.5a) 

wifh C being any eomplex number. Under fhis symmefry, fhe real quantifies 

= and F 5 = |5|2 - (2.5b) 

are invarianf and ean, fhus, be used in fhe eonsfruefion of fhe Kahler pofenfial. The lasf term in fhe 
righf-hand side of fhe equafion for Fs wifh /cs ~ 1 is ineluded in order fo ensure fhaf fhe mass squared 
of S during HI is large and positive - see See. 3.2. If one eombines fhe superpofenfial in Eq. (2.1) 
wifh a eanonieal-like [13] or a logarifhmie Kahler pofenfial involving F_ and Fs, one ean show fhaf 
HI driven by fhe simplesf quarfie pofenfial wifh fransplanekian values of p ean be attained. Namely, in 
Ref. [13], a symmefry similar fo fhe one in Eq. (2.5a) is eonvenienfly applied in fhe ease of fhe SU (2)l 
doublefs of MSSM. This model, fhough, is by now ruled ouf [2] due fo fhe relafively low value of n-s 
(n-g ~ 0.947) and fhe high value of r (r ~ 0.28) fhaf if prediefs - see e.g. Refs. [2,32]. 

We are foreed, fherefore, fo allow a tiny violation of fhe shiff symmefry in Eq. (2.5a) ineluding af 
fhe level of fhe quadrafie ferms in fhe Kahler pofenfial a subdominanf ferm of fhe form 

F+ = |$ + $*|^ (2.6) 

whieh remains invarianf under fhe fransformafion 

4> + C, 5^ S (2.7) 
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coinciding with the one in Eq. (2.5a) only for (7 = 0. Both F_ and F+ respect the symmetries imposed 
on W and generate kinetic mixing between <f> and $ with non-vanishing eigenvalues. However, positive 
eigenvalues of the kinetic mixing of <1> and $ for a logarithmic Kahler potential are provided by F_, 
which has, thus, to play a prominent role - see Sec. 3.1. The quantity F_, contrary to F+, vanishes 
along the trajectory in Eq. (2.4) and it is expected to contribute only to the normalization of the inliaton 
field, whereas F+ is expected to have an impact on both the normalization of the inliaton and the 
inflationary potential. 

Including the terms in Eqs. (2.5b) and (2.6) in a canonical-like Kahler potential, we obtain a model 
which can become just marginally compatible with the data as mentioned in Ref. [13]. Here, we adopt 
two other alternatives, i.e. a purely logarithmic Kahler potential 

A', = -JV. In (l - ^F. + c^F, - ±Fs + ^F^ + . (2.8) 

or a ‘hybrid’ Kahler potential 

K 2 = —N 2 In ^1 — -I- Fs (2.9) 

with one logarithmic term for and $ and one canonical-like kinetic term for S. In both cases, 
positivity of the kinetic energy requires Ni > 0 and N 2 > 0. We also introduced two dimensionless 
coupling constants c_ and c+ with a clear hierarchy c_ S> c+ so that the shift symmetry in Eq. (2.5a) 
is the dominant symmetry of the Kahler potential compared to that in Eq. (2.7). It is worth mentioning 
that our models are completely natural in the ’t Hooft sense because, in the limit c+ —)■ 0 and A —)■ 0, 
the shift symmetry in Eq. (2.5a) becomes exact and a ?7(1) symmetry under which S —)■ (a is a 

real number) appears. 

Note, finally, that the scenario of non-minimal HI investigated in Refs. [17,18] can be recovered 
by doing the substitutions 

c- = 1 + Ncti and c+ = c-ji (2.10) 

with = iVi = 3 in Eq. (2.8) or = A ^2 = 2 in Eq. (2.9) - see below. The symmetries of our 
model, though, prohibit the existence in the Kahler potential of the terms \S\‘^ -|- 

which, generally, violate [32] D-flatness. In this case, the restoration of the D-flatness would require 
the equality of ks^ and kg^, which signals an ugly tuning of the parameters. 

2.3 The Frame Function 

The interpretation of the Kahler potentials in Eqs. (2.8) and (2.9) can be given in the ‘physical’ 
Jordan frame (JF). To this end, we derive the JE action for the scalar fields 2 ;“ = S, <h, We start with 
the corresponding Einstein frame (EF) action within SUGRA [17,32], which can be written as 

S = j d^xs/^ {-]^n + , (2.11a) 

where g is the determinant of the EE metric g^,^, TZ is the EE Ricci scalar curvature, is the gauge 
covariant derivative, and V is the (tree-level) EE SUGRA scalar potential given by 

U = Uf + Ud with Up = e^ (a:“%„f;-- 3|1U|2) and Up = ^< 7 " X]^^.llb) 
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Here, a trivial gauge kinetie funetion is adopted and the summation is applied over the generators Tq of 
the eonsidered gauge group. Also, we use the shorthand notation 




and D, = {Tat Kp with /f, = 


^a/3 ^* 

(2. lie) 


If we perform a eonformal transformation [17,32] defining the JF metrie through the relation 


9fiiy 




we obtain 


9 ^" = -f 

and iz =-^{^Z-a\nn + 3g^^’'^^n^M/2n‘^) . 


(2.12a) 


Here fl is the frame funetion, g is the determinant of TZ is the JF Rieei sealar eurvature and 
Af is a dimensionless parameter whieh quantifies the deviation from the standard set-up [17]. Upon 
substitution of Eq. (2.12a) into Eq. (2.1 la), we end up with the following aetion in the JE 


S = J d‘xV^ Ar + - v'j 

If, in addition, we eonneet QtoK through the following relation 

-n/N = e-^/^ => K = -Nln{-n/N) 


^ 

With V=j;^V. 

(2.12b) 


(2.13a) 


and take into aeeount the definition [17] of the purely bosonie part of the on-shell value of the auxiliary 
field 

A^, = i{Ka,D^z^ - K^D^z*^)/Q, (2.13b) 

we arrive af fhe following aefion 

S = I D.z^D^z*^ - - u) , (2.13e) 

where A^ in Eq. (2.13b) lakes fhe form 


A^ = -iN {^a^D^z^ - n^D^z*^) /6U (2.13d) 

and fhe shorlhand nolalion = U and Uq = U is used. Erom Eq. (2.13a), we ean find fhe 
eorresponding frame funelion during HI as follows 


n _ / (1 + c+(/)2)^i/^ for K = Ki, 

N Eq. (2.4) “1(1 + for K = K 2 , 


(2.14) 


where we look info aeeounl fhe fael lhal F_ = = 0 along fhe direelion in Eq. (2.4). Eqs. (2.13e) 

and (2.14) reveal lhal /t^ represenls fhe non-minimal eoupling lo gravily. Nofe lhal Ihis funelion is 
independenl of c_, whieh is lo be large for HI wilh cj) < I - see below. Seleeling 


N = Ni or N = N 2 for K = Ki or K = K 2 


(2.15) 


respeelively, we ean oblain fhe slandard quadralie non-minimal eoupling funelion. As for Ihe eonven- 
lional ease [17,18] wilh one logarilhm and N = 3, when Ihe dynamies of Ihe fields z" is dominated 
only by the real moduli \z°^\ or when z" = 0 for a / 1 [17], we obtain = 0 in Eq. (2.13d). 
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The only difference w.r.t. the aforementioned conventional case is that now the scalar fields z" have 
non-canonical kinetic terms in the JF due to the term proportional to / This fact does not 

cause any problem, since the canonical normalization of the inflaton retains its strong dependence on 
c_ through 17, whereas the non-inflaton fields become heavy enough during inflation and so they do 
not affect the dynamics - see Sec. 3.1. Furthermore, for M <C mp, the conventional Einstein gravity 
at the SUSY vacuum - in Eq. (2.2b) - is recovered since 

-(U)/Ar~l. (2.16) 

Given that the analysis of inflation in both the EE and JE yields equivalent results [19], we carry out the 
derivation of the inflationary observables exclusively in the EE - see Secs. 3.1 and 3.2. 


3. The Inflationary Set-up 

In this section, we outline the salient features of our inflationary scenario (Sec. 3.1) and then 
present the one-loop corrected inflationary potential in Sec. 3.2. 


3.1 The Tree-level Inflationary Potential 

The linearity of W w.r.t. S allows us to isolate easily the non-vanishing contribution of the infla¬ 
ton to Vp on the inflationary path and avoid the runaway behavior which may be caused by the term 
—3|VFp expiT. Indeed, inserting Eqs. (2.1) and (2.8) or (2.9) into Eq. (2.11b), we find that the only 
surviving contribution to V on the path in Eq. (2.4) is 


Uhio = 


= 


A 2(</>2 - M 2)2 
16 


for K = Ki 


f—N2 

(n 


for K = K 2 . 


(3.1) 


Here we took into account Eqs. (2.14) and (2.15) and the fact that and = f-ji or 

= 1 for K = Ki or K 2 respectively - note that Ksz^^ = 0 for both cases and a = 2 or 3. 
Introducing a new variable n related to the exponents of /t^ in Eq. (3.1), we can cast Uhio in the same 
form for both the K’s in Eqs. (2.8) and (2.9). Indeed, Uhio can be rewritten as 


liiio 


A 2((/)2 - M 2)2 

lg^^(l+n) 


where 


iVi-l = 2(l + n) iK = K, 

N 2 = 2{l + n) \k = K 2 . 


(3.2) 


As anticipated below Eq. (2.7), Uhio depends exclusively on c+ (and not on c_). Given that, during 
HI, (/) S> M and c+(^2 > 1 _ see below -, Uhio and the corresponding Hubble parameter Hui take the 
form 


A20-4»^ 

^HIO - 2(l+n) 

16c^ 




1/2 


HIO 


X(j) 


—2n 


\/3 4V3< 


,1-l-n 


(3.3) 


As a consequence, we obtain an inflationary plateau for n = 0 or a bounded from below chaotic-type 
inflationary potential forn < 0 with cf) in Eq. (2.3) being a natural inflaton candidate. Note that, thanks 
to the shift symmetry in Eq. (2.5a), no mixing term proportional to ks- arises in Uhio in sharp contrast 
with the models of Refs. [8,9], where a similar term (oc ks^) plays a crucial role in achieving large 
values of r in a manner compatible with all observations - see Sec. 4.1. 
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To specify further our inflationary scenario, we have to determine the EF canonically normalized 
fields involved. Note that, along the configuration in Eq. (2.4), the Kahler metric defined in 
Eq. (2.11c) takes, for both choices of K in Eqs. (2.8) and (2.9), the form 


1 


i^ap) = diag (Mii-,iT5s*), where Mk = 

Jn 


K K 


K K 


with K = c-Zti — iVc+, K = Nc 


,2 J.2 


(3.4) 

and N defined in Eq. (2.15). Given that Kss* = or 1 for iT = Ki or K 2 respectively, the 
canonically normalized components 's,s of S - see Eq. (2.3) - are defined as follows: 


(S,S) = y/KsS*{s,s). 


(3.5) 


The matrix Mk can be diagonalized via a similarity transformation involving an orthogonal matrix Uk 
as follows: 


UkMkUk = diag (k+, «:_), where Uk = 


1 

71 


1 1 
-1 1 


(3.6) 


and the eigenvalues of Mk are found to be 

c_(l — Nr±) + c+c_(l + Nr±)(fP‘ 


K+ = 


c_ , C-{1—Nr±) 

~ —— and K- = - 


c_ 


fl fn fn fn 


- (3.7) 


where the approximate results hold for <C f/N and positivity of k- can be assured only if 

r± < 1/N with = c+/c_, 


(3.8) 


i.e. if c_ > Nc+, as we anticipated below Eq. (2.7). This fact has to be contrasted with the original 
scenario of non-minimal HI [17,18], where such a constraint is not necessary - as can be verified by 
inserting Eq. (2.10) into Eq. (3.7). In our present cases, the kinetic terms for 2;“ = $, can be brought 
into the following form 


K — — 


X2 


1 




+ 


/I . 




1 


■2 


+ 0_i_ + 6_ + 9 


$ ) 


(3.9a) 


where 9± = ((9 ± 9) /\/2 and the dot denotes derivation w.r.t. the cosmic time t. In the last step, 
we introduce the EF canonically normalized fields, which are denoted by hat and can be obtained as 
follows: 


d(p 


— J — y77+, 


J(t)9+ 

~72 


, e_ = 





9_, and 9q, = ( 9q, — 


TT 


(3.9b) 


Note, in passing, that the spinors 'il)s and associated with the superfields S and <1>± = (4> ± 4>)/\/2 
are normalized similarly, i.e. ips = '\/Kss*'4’s and Integrating the first equation in 

Eq. (3.9b), we can express the canonically normalized EF real field (/> as follows: 


— <Pc + 




arcsinhyTi+i;/* 


(3.10) 



with (f)c being a constant of integration, which we take equal to zero. Note that (p is practically inde¬ 
pendent of N (and n) - see Eq. (3.7). Solving Eq. (3.10) w.r.t. 4>, we can express Fhio in Eq. (3.3) in 
terms of (p as follows: 


f^HIO — 


tanh^ ^J^(p 
4c+ cosh^” y/f±(p 


(3.11) 


Eor n = 0 and A = c+ = 1, Vhio coincides with the potential encountered in the so-called T-models [6] 
arising from the spontaneous breaking of (super)conformal invariance. We observe that, although /t^ 
in Eq. (2.14) and Vhio in Eq. (3.11) are independent of c_, (p depends heavily on c_ and, therefore, 
it can be much larger than (p facilitating the attainment of HI with subplanckian values of (p. As a 
consequence, the initial fields <h and $ - see Eq. (2.3) -, which are closely related with cp, can remain 
also subplanckian, as required for a meaningful approach to SUGRA. 


3.2 Stability and One-loop Radiative Corrections 

To ensure the validity of our inflationary proposal, we have to check the stability of the direction 
in Eq. (2.4) w.r.t. the fluctuations of the fields which are orfhogonal fo Ibis direction, i.e. we have fo 
examine fhe fulfillmenl of fhe following condifions: 


dV 


dx^ 


= 0 and fnia > Q wifh s, s. 


(3.12a) 


Eq. (2.4) 


Here fri^a are fhe eigenvalues of fhe mass-squared mafrix wifh elemenfs 



d'^V 

dx^^dx^ 


wifh 

Eq. (2.4) 


0_, 0$, s, s. 


(3.12b) 


Diagonalizing we consfrucf fhe scalar mass-squared specfrum of fhe fheory along fhe direcfion 
in Eq. (2.4). In Table 2, we presenf approximate expressions of fhe relevanf masses squared, which 
are quite close fo fhe rafher lengfhy exacf expressions. Note, however, fhaf our numerical compulation 
uses fhe exacf expressions. In Ihis fable, we also include fhe mass squared of cp as well as fhe 
masses squared of fhe chiral fermions, fhe gauge boson Abl, and fhe gaugino Xbl which are used in 
our analysis below. 

Erom fhe formulas displayed in Table 2, we can infer fhaf fhe sfabilify of fhe pafh in Eq. (2.4) is 
assured since Eq. (3.12a) is fulfilled. In particular, if is evidenf fhaf > 1 assisfs us fo achieve > 0 
for K = Ki - in accordance wifh fhe resulfs for similar models in Refs. [8,17]. On fhe ofher hand, 
for K = K2, ml > 0 even wifh ks = 0. However, since fhere is no observafional bin! [2] for large 
non-Gaussianify in fhe cosmic microwave background, we should make sure fhaf all fhe ’s for fhe 
scalar fields in Table 2 excepf are greafer fhan 77yj during fhe lasf 50 — 60 e-foldings of HI. This 
guaranfees fhaf fhe observed curvafure perfurbafion is generafed wholly by (p as assumed in Eq. (4.4) 
- see below. Requiring fhaf ml >> enfails fhe existence of a non-vanishing ks for K = K2 
too. Due fo fhe large effective masses fhaf fhe scalars acquire during HI, fhey enfer a phase of damped 
oscillafions abouf zero. As a consequence, fhe (p dependence in fheir normalization - see Eq. (3.9b) - 
does nol affecl fheir dynamics. 
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Fields 

Eigenstates 

Masses Squared 




K = Ki 

K = K2 

3 real scalars 


m2 

^^Hl 


0+ 

m^+ 

6(1 - l/Ni)Hl, 

ei/^i 


0<j) 


^BL + 6 (^1 - 7 ^^ 77hi 

^BL + 

1 complex scalar 

s, s 

— 2 
m. 

Q{2ksfn-^)Hl, 

^‘2ksHyii 

1 gauge boson 

^BL 

^BL 

p2c_(l - Nr±)(t?/f'R, 

4 Weyl spinors 



6(2+c+(3—u2 

^Hl 

6(2+c+( 2-A2)</.2)2 f>2 

"hi 



^BL 

p2c_(l - Nr±)(t)^/fn 


Table 2: The mass-squared spectrum for K = Ki and K = K 2 along the inflationary trajectory in Eq. (2.4) for (j> 1. 

N is defined in Eq. (2.15) and fj is given by Eq. (4.10) - see below. To avoid very lengthy formulas, we neglect terms 
proportional to M <C 0. 

Considering SUGRA as an effective theory below mp allows us to use the well-known Coleman- 
Weinberg formula [23] in order to find the one-loop corrected inflationary potential 

Vui = Vhio + AVhi with AVri = ’ (3.13a) 

i 

where the sum extends over all helicity states i of the fields listed in Table 2, Fj is the fermion number 
and inf the mass squared of the ith helicity state, and A is a renormalization mass scale. The consistent 
application of this formula requires that the rhf’s which enter into the sum are: 

• Positive. As a consequence and following the common practice [8,9,20], we neglect the con¬ 
tribution of to A1 /hi since this mass squared turns out to be negative in the largest part of 
the parameter space of our models. Let us recall here that Eq. (3.13a) is valid only for a static 
configuration and is uniquely defined only in the extremum points. The proper calculation should 
use the time-dependent background - see e.g. Ref. [33]. For the non-minimal inflation, this is 
not done up to now, but we do not expect to get a very unexpected effect if the full computation 
is consistently carried out. 

• Much lighter than a momentum cut-off squared, which is here considered as large as mp. As a 
consequence, we do not take into account the contributions from and inf^ to APhi since 
these masses squared are much larger than mp in our case given that c_ S> 1 - see below. This 
stems from the form of k_ in Eq. (3.7) and does not occur in the case of the standard non-minimal 
HI [18,20], as can be checked by substituting Eq. (2.10) into the expressions for these masses 
squared in Table 2. 

Having in mind the above subtleties and neglecting contributions from the gravitational sector of the 


10 




theory, the one-loop RCs read 


AFhi = 


+ 2m^ln^-4m^^ In ^ 


(3.13b) 


The renormalization seale A ean be determined by requiring [24] that = 0 or ) = 0. 

Let us, finally, stress here that the non-vanishing value of (p during HI breaks spontaneously 
U{1)b-l leading to a Goldstone boson 6-. This is ‘eaten’ by the gauge boson Abl, whieh then 
beeomes massive. As a eonsequenee, six degrees of freedom before the spontaneous breaking (four 
eorresponding to the two eomplex sealars <I> and $ and two eorresponding to the massless gauge boson 
Abl U{1)b-l) are redistributed as follows: three degrees of freedom are assoeiated with the real 
propagating sealars {p, and 9^), whereas the residual one degree of freedom eombines together 
with the two ones of the initially massless gauge boson Abl^o make it massive. From Table 2, we ean 
deduee that the numbers of bosonie (eight) and fermionic (eight) degrees of freedom are equal, as they 
should. 


4. Constraining the Parameters of the Models 

We will now outline the predictions of our inflationary scenarios in Sec. 4.2 and confront them 
with a number of criteria introduced in Sec. 4.1. 

4.1 Inflationary Observables - Constraints 

Our inflationary settings can be characterized as successful if they can be compatible with a number 
of observational and theoretical requirements, which are enumerated in the following - cf. Ref. [34]: 

4.1.1. The number of e-foldings 



that the pivot scale fc* = 0.05/Mpc suffers during HI has to take a certain value to resolve the horizon 
and flatness problems of standard hot big bang cosmology. This requires [2] that 


A* ~ 61.5-1- In 


1 - Srui-h T^‘^grh*TX 

12(1-|-rUi-h) 30VHi(i;Af) 


In — ^ , - — Infifrh*, 


where we assumed that HI is followed, in turn, by a phase of damped inflaton oscillations with mean 
equation-of-state parameter rurh, a radiation dominated era, and a matter dominated period. Here, Tj-h 
is the reheat temperature after HI, ^(rh* is the energy-density effective number of degrees of freedom at 
Trh - for the MSSM spectrum, we take ^(rh* = 228.75 -, 0* [(/>*] is the value of p [p] when /c* crosses 
outside the inflationary horizon, and pf [pf\ is the value of p [p] at the end of HI, which can be found, 
in the slow-roll approximation, from the condition 


max{e((/>f), \g{pi)\} = 1 


(4.3a) 
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with the slow-roll parameters ealeulated as follows: 


1 / _ 1 / \ 
2 V ^Hi y ~ 2J2 J 


and 




1 j Hit,(/)(/) 

^ V ^HI 




(4.3b) 


Given that, for a power-law potential (j)"', we have [34, 35] Wj^h = {n — 2)/{n + 2), we take for our 
numeries Wrh = 1/3, whieh eorresponds preeisely to n = 4. Although we expeet that, in the our eases, 
rurh will deviate slightly from this value, we eonsider this value quite reliable sinee, for low values of 
(j), our inflationary potentials ean be well approximated by a quartie potential. As a eonsequenee, our 
set-up is largely independent from Tj-h - see Eq. (4.2). 


4.1.2. The amplitude As of the power speetrum of the eurvature perturbation generated by cp at the 
pivot seale /c* must be eonsistent with the data [2]: 


A//2 = 


1 

IV#*)I 2v/37r \VniA<P*)\ 


4.627 X 10 


-5 


(4.4) 


where we assume that no other eontributions to the observed eurvature perturbation exist. 


4.1.3. The remaining inflationary observables, i.e. the sealar speetral index rig, its running Og, and the 
tensor-to-sealar ratio r, whieh are given by 


(a) ng = 1 - 6e* + 2r)*, (b) Og 



(ns - 1)^) - 2^*, and (c) r = 16e*, (4.5) 


where 




T/2 


I4ti,0 v,<i> 

V/n 


+ 2rie 


(4.6) 


and the variables with subseript * are evaluated ni(l) = (p*, must be in agreement with the fitting of the 
data [2] with the ACDM-|-r model, i.e. 


(a) Us = 0.968 ± 0.009 and (b) r < 0.12, (4.7) 

at 95% e.l. with |ag| <C 0.01. Although eompatible with Eq. (4.7b), the present eombined Planck and 
^iCEvHKeck Array results [1] seem to favor models with values of r of order 0.01 - see Eq. (1.1). 

4.1.4. To avoid eorreetions from quantum gravity and any destabilization of our inflationary seenario 
due to higher order non-renormalizable terms, we impose two additional theoretieal eonstraints on our 
models - keeping in mind that Vm{(pi) < Vm{(p-iP)\ 

(a) < 1 and (b) (/* < 1. (4.8) 

As we will show in See. 5, the UV eutoff of our model is equal to unity (in units of mp) for r± < 1 
and so no eoneerns regarding the validity of the effeetive theory arise. 
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4.1.5 The U (1) b-l gauge symmetry does not generate any extra eontribution to the renormalization 
group running of the MSSM gauge eoupling eonstants and so the seale M and the relevant gauge 
eoupling eonstant qb-l can be mueh lower than the values dietated by the uniheation of the gauge 
eoupling eonstants within the MSSM. However, for definiteness, we eonsider here the most predietive 
ease in whieh qb-l = 9 id ~ 0.7 is the SUSY GUT gauge eoupling eonstant) and M is determined by 
requiring that (<h) and ($) take values eompatible with the uniheation of the MSSM gauge eoupling 
eonstants. In partieular, the SUSY GUT seale Mgut — (2/2.433) x 10“^ is to be identihed with the 
lowest mass seale of the model at the SUSY vaeuum in Eq. (2.2b), i.e. 


Vc (1 - Nr±)g M ^ ^ M = 

\l9‘^c-{l - Nr±) - 


(4.9) 


The requirement that the expression g^c-{l — Nr±) — c+Mq^t is positive sets an upper bound on c+ 
for every c_. Namely, we should have c+ < 5^(1 — iVr±)c_/MQUT’ whieh, however, is too loose to 
restriet the parameters. 


4.2 Analytic Results 


Our analytie results are based on the tree-level inhationary potential in Eq. (3.2) and are identieal 
for both K = Ki and K = K 2 provided that Ni and N 2 are related to n as shown in this equation. 
Note that, not only the form of Uhio in Eq. (3.2), but also the eanonieal normalization of (p is praetieally 
identieal in the two eases - see Eqs. (3.7) and (3.9b). The slow-roll parameters read 


8(1 — nc+(/)^)^ 


and rj = A 


3 + c+(/>^ (n (Anc+cfP' — 9) — 2) 

c-P^fn 


(4.10) 


The termination of HI is triggered by the violation of the rj eriterion at a value of p equal to p{. Sinee 
pi <C (/>*, the slow-roll parameters in Eq. (4.10) ean be well approximated by performing an expansion 
for small values of p. We hnd 


, 1 — (1 + 2n)c+0^ 


C-C 


and ~ 4 


3 - (5 -f 9n)c+p‘^ 
C-p'^ 


Employing these expressions, p{ is ealeulated to be 


rj (Pi) = 1 ^ (/)f ~ 2-spP ^c_ -f 20c+ -f 36nc+^ 


-1/2 


(4.11) 


(4.12a) 


Note that the violation of the ? eriterion oeeurs at p = pf sueh that 


e(0f) 


1 ^ ~ 2 V 2 ^c_ -|- 8c+ + 16nc+^ < p{. 


(4.12b) 


We proeeed with our analysis presenting separately our results for the two radieally different eases: the 
ease n = 0 in See. 4.2.1 and the ease n < 0 in See. 4.2.2. 


4.2.1 The n = 0 Case 

Given that pi <C 0*, iV* ean be ealeulated via Eq. (4.1) as follows: 

A* ~ c_ {Pi - pf) /8 ^ ~ 2Y^2iV*/c_. (4.13a) 
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Obviously, HI with subplanckian values of (p can be attained if 


0* < 1 ^ c_ > 8A^* ~ 480 


(4.13b) 


for A/* ~ 60. Therefore, large values of c_ are dictated, whereas c+ remains totally unconstrained by 
this requirement. Replacing Vhio from Eq. (3.2) in Eq. (4.4), we find 


4 1/2 _ 

® 327rY3 + 3c+02 



(4.14) 


where /o* = fn{(p*) = 1 + (for n = 0) and Eq. (4.13a) was use in the last step. Inserting, 

finally, fhis equation info Eq. (4.5), we find fhe following expressions for ns, Os, and r: 

2 1 1 - 2/o* - 2fl 16 

ns 1 —^ -, Os ^and r ~ -• (4.15) 

N^fo* iV*/o* iV*/o* 

Therefore, a clear dependence of fhe observables on r± arises. If is worth noticing fhaf fhese resulfs 
coincide wifh fhe ones obfained for fhe model of kinefically modified non-minimal inflation esfablished 
in Ref. [10] wifh m = 0 and n = 4 - in fhe nofafion of fhis reference. 


4.2.2 The n < 0 Case 

When n < 0, W* can be esfimafed again fhrough Eq. (4.1) wifh fhe resulf 


1 

8nr± 1 — ’ 


(4.16) 


which, obviously, cannof be reduced fo fhe one found for n = 0 - cf. Eq. (4.13a). Neglecting pi in fhe 
lasf equably - since (/>f <C </>* - and solving w.r.f. (/>*, we find 


cii , -—where Cn = e . (4.17a) 

V 

Nofe fhaf fhe dependence of p^ on W* is radically differenf from fhe one in Eq. (4.13a) and resembles 
fhe dependence found in Refs. [8,9] for n < 0. However, (/>* can again fulfill Eq. (4.8b) since 

p^<l => c_ > (4.17b) 

\n\r± 

where fhe lower bound on c_ furns ouf lo be r±-dependenl - in confrasl fo fhe case of Eq. (4.13b). 
Subslifuling Eq. (4.17a) info Eq. (4.4) and solving w.r.f. A, we end up wifh 

A ~ 32TTy/^senfn*^^^c_{nr±/{l - e„))^/^, (4.18) 


where /„* = = (1 + n — en)ln for n < 0. We remark fhaf A remains proportional lo c_ (for 

fixed n and r±) as in fhe case wifh n = 0 - cf. Eq. (4.14). Inserting Eq. (4.17a) info Eq. (4.10) and 
employing Eq. (4.5a), we find 


8r± (2ne^ + 2(1 + n) - (2 + n)en) 

ns - 1 +-^- 7z -TT-^ 

(1 (^n)jn* 


(4.19a) 


14 



From this expression, we see that n < 0 and r± <1 assist us to reduee n-s so as to beeome eonsiderahly 
lower than unity as required by Eq. (4.7a). Using Eqs. (4.17a), (4.10), and (4.5b, c), we arrive at 

64enr^ (e^(2 + n)(2n - 1) - (1 +n)(2 + n) - 4en(n^ - 1)) , 128ne^r± 

=-n-1272- ^ - n - \T~ ' 

(1 ^n) /n* (1 

(4.19b) 

Erom the last result, we eonelude that mainly the faet that |n| / 0 and seeondarily the faet that n < 0 
help us to inerease r. 

4.3 Numerical Results 

Adopting the definition of n in Eq. (3.2), our models, whieh are based on W in Eq. (2.1) and K in 
Eq. (2.8) or (2.9), ean be universally deseribed by the following parameters: 

A, n, c_, c+, ks, /c_, and ks-- 

Reeall that A* turns out to be independent of T^h, as explained in See. 4.1, and M, whieh is determined 
by Eq. (4.9), does not affeet the inflationary dynamies and the predietions sinee M <C </> during infla¬ 
tion. Erom the remaining parameters, ks influenees only in Table 2 and ks- enters only into the 
higher order terms - not shown in the formulas of Table 2 - in the expansions of rh‘^_^ and fhg^. On 
the other hand, k- does not appear at all in our results. Given that the eontribution of AUhi to Uhi in 
Eq. (3.13b) ean be easily tamed with a suitable seleetion of A - see Table 3 below -, our inflationary 
outputs are essentially independent of these three parameters, provided that we ehoose them so that the 
relevant masses squared are positive. To ensure this, we set ks = ks- = 1 throughout our ealeula- 
tion. Moreover, the bulk of our results are independent of the ehoiee between K = Ki or K = K 2 , 
espeeially for r± <0.1 sinee J in Eq. (3.9b) remains undistinguishable. However, for definiteness, we 
present our results for K = Ki, unless otherwise stated. 

Eor fixed n, the remaining three free parameters of our models during HI, whieh are c_, c+, and 
A, ean be redueed by one leaving us with the two free parameters r± and A/c_. This faet ean be 
understood by the following observation: If we perform the rescalings 

and 5^5, (4.20) 

the superpotential W in Eq. (2.1) depends on A/c_ (M is not important as we explained) and the Kahler 
potential K in Eq. (2.8) or (2.9), for 5=0, = $*, and fixed n, depends on r±. As a consequence, 

Uhio depends exclusively on A/c_ and r± via f-ji in Eq. (2.14). The confrontation of these parameters 
with observations is implemented as follows: Substituting Uhi from Eq. (3.13a) in Eqs. (4.1), (4.3b), 
and (4.4), we extract the inflationary observables as functions of n, r±, A/c_, and (/>*. The two latter 
parameters can be determined by enforcing the fulfillment of Eqs. (4.2) and (4.4), whereas n and r± 
largely affect the predictions for rig and r and are constrained by Eq. (4.7). Moreover, Eq. (4.8b) bounds 
c_ from below, as seen from Eqs. (4.13b) and (4.17b). Einally, Eq. (3.8) provides an upper bound on 
r±, which is different for N = Ni and N = N 2 , discriminating slightly the two cases. 

We start the presentation of our results by checking the impact of AUhi in Eq. (3.13b) on our 
inflationary predictions. This is illustrated in Table 3, where we arrange input and output parameters 
of our model with K = Ki which are consistent with the requirements of Sec. 4.1. Namely, we fix 
r± = 0.015 and n = —1/50 orn = —1/20, which are representative values as seen from Eig. 1 below. 
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Input Parameters 

— 1/n 

50 

20 

50 

20 

50 

20 

C-/IO2 

5.03 

22.3 

5.03 

22.3 

5.03 

22.3 

r± 

0.015 

0.015 

0.015 


1 

0.5 

1 

0.5 

1 

0.5 

Output Parameters 


0 

II 

<1 

AFhi(0 

*) = o 

AVHi((('f) = 0 

A/10-3 

1.13 

5.21 

1.13 

5.21 

1.13 

5.21 

(/.f/O.l 

1.4 

0.67 

1.4 

0.67 

1.4 

0.67 

A/10-3 

- 

- 

20.1 

23.5 

1.02 

1.07 

A* 

58.5 

58.1 

58.5 

58.1 

58.5 

58.1 

ris/O.l 

9.66 

9.68 

9.66 

9.68 

9.66 

9.68 

1 

0 

1—1 

1 

6.3 

6.4 

6.3 

6.4 

6.3 

6.4 

r/0.01 

3.8 

4.77 

3.8 

4.77 

3.8 

4.77 


Table 3: Input and output parameters of the model which are compatible with all the requirements of Sec. 4.1 for 
ks = ks- = 1- We use the tree-level potential by switching off the RCs (i.e. taking AVhi = 0) or the one-loop 
corrected potential with the renormalization scale A determined such that = 0 or AVHi(</>f) = 0, as 

indicated in the table. 

In the second and third columns of this table, we accumulate the predictions of the model with the 
RCs switched off, whereas, in the next columns, we include AChi- Following the strategy adopted in 
Ref. [24], we determine A by requiring AChi(</>*) = 0 or AyHi(';^f) = 0. We can easily deduce that our 
results do not change after including the RCs with either determination of A, since AFhi remains well 
suppressed in both cases. Note that the resulting A is well below unity in the two cases with its value 
in the case with AFhi(</>*) = 0 being larger. Therefore, our findings can be accurately reproduced by 
using Fhio instead of Fhi- This behavior persists even if we take K = K 2 . In this case, AFhi assumes 
even lower values, especially if we select A such that AkHiCi^*) = 0. This is due to the fact that the 
values of m* entering into Eq. (3.13b) for K = K 2 different from those for K = Ki. However, 
this does not cause any differentiation between the two models. 

The predictions of our models can be encoded in the ns — ro .002 plane, where ro .002 is the value 
of r at the scale k = 0.002/Mpc. This is shown in Fig. 1 for n = 0 (solid line), n = —1/50 
(dashed line), and n = —1/20 (dot-dashed line). The variation of r± on each line is also depicted. 
To obtain an accurate comparison with the marginalized joint 68% [95%] regions from the Planck, 
Bicep2I Keck Array and BAO data - they are also depicted as dark [light] shaded areas -, we compute 
^ 0.002 = 16e'((^0.002)> where 0o.oo2 is the value of cj) when the scale k = 0.002/Mpc, which undergoes 
-^ 0.002 = A* + 3.22 e-foldings during HI, crosses outside the horizon of HI. For n = 0 and r± < 1, the 
line is almost straight and essentially coincides with the corresponding results of Refs. [10,22]. For low 
values of r±, this line converges toward the values of ns and ro .002 obtained within the simplest model 
with a quartic potential, whereas, for larger values of r±, it crosses the observationally allowed corridors 
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Figure 1: Predicted curves in the ris — ro,oo 2 plune torn = 0 (solid line), n = —1/50 (dashed line), n = —1/20 
(dot-dashed line), ks = ks- = 1, and various values of r± indicated on the curves. The marginalized joint 68% 
[95%] regions from the Planck, BlCEP2/Keck Array and BAO data are depicted as dark [light] shaded areas. 

approaching its universal attractor value [22] for 3 > 1. We cut this line at ~ 1/3 [r± ~ 1/2] 
for K = Ki [K = K 2 ], where the hound in Eq. (3.8) is saturated. This hound overshadows the one 
derived from the unitarity constraint - see Sec. 5 - and restricts ro .002 to be larger than 0.0028 [0.0019] 
for K = Ki [K = K 2 ]. For quite small values of r±, the curves corresponding to n < 0 converge to 
the curve for n = 0. However, for larger values of r±, these curves move away from the n = 0 line 
turning to the right and spanning the observationally allowed ranges with quite natural values of r±, 
consistently with Eqs. (4.19a) and (4.19b), which are in excellent agreement with the numerical results. 
Similarly to the n = 0 case, the n < 0 cases too provide us with a lower bound on r. Specifically, for 
n = —1/50 [n = —1/20], we obtain ro .002 ^ 0.0123 [ro.oo 2 > 0.03]. Finally, we remark that, for any 
n, we can define a minimal r™™ and a maximal value of r± in fhe marginalized join! 95% region. 
Specifically, we find 


^min _ 4 5 ^ ;^q-3 

and 

^max _ 4/3 [ 4 / 2 ] 

for 

n = 

0 , 


^min _ 4_g ^ ;^q-3 

and 

^max _ Q 43 

for 

n = 

-1/50, 

(4.21) 

~ 5.0 X 10-3 

and 

^max _ Q Q52 

for 

n = 

- 1 / 20 , 



where fhe values given for n = 0 correspond to K = Ki [K = K 2 \. Nofe fhaf increasing \n\, 
decreases and becomes more nafural according fo fhe argumenf below Eq. (2.9). 

The sfrucfure of Ihi as a function of cj) for 1 /)* = l,r± = 0.015, and fhe values of n employed in 
Fig. 1 is displayed in Fig. 2. The values of n, A, and r±, shown in fhis figure, yield rig = 0.964,0.966, 
or 0.968 and r = 0.033,0.038, or 0.047 for n = 0, —1/50, or —1/20 respectively. The corresponding 
values of c_ are (4.66, 5.03, or 5.6) x 10^, whereas fhe values of 0* derived from Eq. (3.10) are 13.87, 
14.17, or 14.6. We verify fhaf observable values of r are associated wifh fransplanckian values of 
1 /)* in accordance wifh fhe Fyfh bound [36]. This fad, fhough, does nol invalidate our scenario since 
fhe corresponding values of fhe inifial inflalon cf), which is direcfly relafed fo fhe superfields 4> and 4> 
appearing in fhe definifion of our models in Eqs. (2.1) and (2.8) or (2.9), remain subplanckian - cf. 
Eq. (4.8b). We also remark fhaf, in all cases, %j]j^(i/)*) furns ouf fo be close fo fhe SUSY GUT scale 
Mgut — 8.2 X 10“^, which is imperative - see e.g. Ref. [37] - for achieving values of r close fo 0.1. 
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Figure 2: The inflationary potential Vhi as a function of (p for (^ > 0, r± ~ 0.015, and n = 0, A = 1.04 x 10“^ 
(light gray line), n = —1/50, A = 1.14 x 10“^ (black line), or n = —1/20, A = 1.3 x 10“^ (gray line). The 
values of and pi are also indicated. 




X 


X 


(a) 


(b) 


Figure 3: Allowed regions in the A — c_ (a) and A — Og (b) plane for ks = ks- = 1 and n = 0, —1/50, or 
—1/20 bounded by light, dark, or normal gray lines, respectively. The conventions adopted for the various lines 
are shown in the panel (a). 


We finally observe that the slope of Vm elose to (/> = (/>* inereases with \n\. This is expeeted to elevate 
?- see See. 4.2 - and, via Eq. (4.5c), r too. 

To speeify further the allowed ranges of the parameters of our models, we plot in Figs. 3-(a) and 
-(b) the allowed regions in the A — c_ and A — Og planes, respeetively. The eonventions adopted for 
the various lines are shown in Fig. 3-(a). In partieular, the boundary eurves of the allowed region for 
n = 0, —1/50, or —1/20 are represented by light, dark, and normal gray lines, respeetively. The 
dot-dashed and dashed lines eorrespond, respeetively, to the minimal and maximal values of r± given 
in Eq. (4.21), whereas, on the thin short-dashed lines, the eonstraint of Eq. (4.8b) is saturated. The 
perturbative bound A < 3.5 (so that the expansion parameter A^/dvr < 1) limits the various regions 
at the opposite end by a thin solid line. For K = K 2 , the light gray dashed line is expeeted to be 
transported to lower values of c_, but keeping the same slope. Note that the dot-dashed lines eoineide 


18 












1 



10 20 30 40 50 60 70 80 90 100 


-1 /n 

Figure 4: Allowed (shaded) region compatible with Eqs. (4.2), (4.4), (4.7), and (4.8) in the (— 1/n) — r± plane 
for ks = ks- = 1- The conventions adopted for the various lines are also shown. 

with each other in Fig. 3-(a) and this is almost the case with the dark and light gray dot-dashed lines in 
Fig. 3-(b) too. We observe that A/c_ remains constant for fixed r±, as expected by the argument below 
Eq. (4.20). The required c_ for subplanckian excursions of tp - overall, we obtain 0.014 <(/)<!- 
is quite large and increases with |n|. On the other hand, |as| remains sufficiently low. Therefore, our 
models are consistent with the fitting of the data with the ACDM-i-r model [2]. 

Concentrating on the most promising cases with n < 0, we delineate, in Fig. 4, the overall allowed 
region of our models by varying continuously n and r±. The conventions adopted for the various lines 
are also shown in the figure. In parficular, fhe dashed [dof-dashed] line corresponds fo Ug = 0.977 
[n-s = 0.959], whereas fhe solid line is obfained by fixing ng = 0.968 - see Eq. (4.7a). On fhe fhin 
solid boundary line, fhe bound in Eq. (4.7b) is safurafed. We remark fhaf, as r± increases wifh fixed n, 
Hg increases too, while r decreases. This agrees wifh our findings in Eig. 1. Nofe fhaf, for n > —1/30, 
fhe fhick dof-dashed and fhe fhin solid lines coincide. Overall, for Ug = 0.968 and ~ 58, we find: 

1.86 < 10®— < 3.6 wifh 8 < -- < 100, 1 < ^ < 3.6 , and 1.8 < < 12 . (4.22) 

Erom fhese resulfs, we infer fhaf r± lakes more nalural - in fhe sense of fhe discussion below Eq. (2.9) 
- values (lower lhan unify) for larger values of |n|. 

5. The Effective Cut-off Scale 

An oufslanding Irademark of our selling is fhaf perlurbalive unilarily is relained up lo mp, despife 
fhe facl fhaf fhe implemenlalion of HI wifh subplanckian values of (p requires relalively large values of 
c_ - see Eqs. (4.13b) and (4.17b). To show Ihis, we exlracl fhe UV cul-off scale Auv of fhe effeclive 
Iheory following fhe syslemalic approach of Ref. [27]. Eel us firsl clarify fhaf, allhough fhe expansions 
presenled below aboul (0) = M <C 1 are nol valid [26] during HI, we consider fhe resulling Auv 
as fhe overall cul-off scale of fhe Iheory, since rehealing is regarded as an unavoidable slage of fhe 
inflafionary dynamics [27]. 
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The canonically normalized inflaton can be written as follows - see Eq. (3.9b): 


5(J) = {J)64i with S(t> = (t>-M and (J) = y/( k+) ~ \/ C-Kfn), (5.1) 

where the last (approximate) equality is valid only for r± <C l/N - see Eq. (3.8). Note, in passing, 
that the mass of dip at the vacuum is calculated to be 

/- \V2 /- „\l/2 am 

^ -^====== (5.2) 

with numerical values (1.5 — 3.3) x 10“® along the solid line of Eig. 4 where ns ^ 0.968. Given 
that (/t^) ~ 1, |n| <C 1, and A/c_ is fixed - see Eqs. (4.14) and (4.18) -, the variation of fhs^ is 
mainly generated by the variation of r±. In other words, for fixed n, fhsip depends only on r± and 
nof on X,M,C-, or c+ separafely. Nofe fhaf fhe resulting values of are almost two orders of 
magnitude lower than the values obtained in similar models [8,9,18] and so the successful activation 
of the mechanism of non-thermal leptogenesis [28] is an important open issue. 

The fact that 5(j) does not coincide with 5(j) in Eq. (5.1) - contrary to the case of standard non- 
minimal HI [25,26] - ensures that our models are valid up to mp = 1 - cf. Ref. [10]. Taking into 
account Eq. (3.9a) and calculating the action S in Eq. (2.11a) on the path of Eq. (2.4) with spatially 
constant values of (p, we find 

S = y \ -) ’ 

where J is given in Eq. (3.9b) and Hhio in Eq. (3.2). Expanding abouf (p 
expressing fhe resulf in terms of cp (~ 5(p) using Eq. (5.1), we obtain 

, 2 

j‘^(p‘^ = — r±^ + — 5Nr^(p^ + • • . (5.3b) 

If we insert Eq. (2.10) into Eqs. (3.7) and (3.9b) and repeat the same expansion, we get the notorious 

term 2Nc^(p‘^ [25,26], which entails Apv ^ 1 for large ct^. In our case, however, cti is replaced 

by r± which remains low enough - although c+ and c_ may be large - rendering the effective theory 
unitarity safe. Expanding similarly Vhio in terms of <p, we have 

f^Hio = (l - 2(1 + n)r±^ + (3 + 5n + -j • (5.3c) 

This expression, for n = 0, reduces to the one presented in Ref. [10], whereas the expression for J in 
Eq. (5.3b) differs slightly from the corresponding one in this reference, due to the different normaliza¬ 
tion of (p in Eqs. (3.7) and (3.9b). Since the positivity of k- in Eq. (3.7) entails r± < l/N < 1, our 
overall conclusion is that our models do not face any problem with perturbative unitarity up to mp. 

6. Conclusions 

We presented models of Higgs inflation in SUGRA, which accommodate inflationary observables 
covering the ‘sweet’ spot of the recent joint analysis of the Bicep 2/Keck Array and Planck data. Our 


(5.3a) 

= {(p) = M ~ 0 and 
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models, at the renormalizable level, are tied to a unique superpotential determined by an R and a gauge 
U{1)b-l symmetry. We seleeted two possible Kahler potentials K, one logarithmie and one semi- 
logarithmie - see Eqs. (2.8) and (2.9) -, whieh respeet the symmetries above and a mildly violated shift 
symmetry. Both K’s lead to praetieally identieal inflationary models. The key-point in our proposal 
is that the eoeffieient c_ of the shift-symmetrie term in the Kahler potentials does not appear in the 
supergravity inflationary sealar potential expressed in terms of the original inflaton field, but strongly 
dominates the eanonieal normalization of this inflaton field. The inflationary seenario depends essen¬ 
tially on three free parameters (n, A/c_, and r±), whieh are eonstrained to natural values leading to 
values of rig and r within their I — a observational margins. Indeed, adjusting these parameters in 
the ranges n = —(0.1 — 0.01), A/c_ = (1.86 — 3.6) x 10“®, and r± = 0.01 — 0.036, we obtain 
ris — 0.968 and 0.018 < r < 0.12 with negligibly small values of |as|. Imposing a lower bound on 
c_, we sueeeeded to realize HI with subplanekian values of the original inflaton, thereby stabilizing 
our predietions against possible higher order eorreetions in the superpotential and/or the Kahler poten¬ 
tials. Moreover, the eorresponding effective theory remains trustable up to mp. We also showed that 
the one-loop RCs remain subdominant for a convenient choice of the renormalization scale. Finally, 
the scale of U{1)b-l breaking can be identified with the SUSY GUT scale and the mass fhs^ of the 
normalized inflaton is confined in the range (1.5 — 3.3) x 10“®. 

As a last remark, we would like to point out that, although we have restricted our discussion to the 
Gb-l gauge group, the HI analyzed in this paper has a much wider applicability. It can be realized 
within other SUSY GUTs too based on a variety of gauge groups - such as the left-right, the Pati- 
Salam, or the flipped SU{5) group - provided that a conjugate pair of Higgs superfields is used in 
order to break the symmetry. In these cases, the inflationary predictions are expected to be quite similar 
to the ones discussed here. The discussion of the stability of the inflationary trajectory may, though, be 
different, since different Higgs superfield representations may be involved in implementing the GUT 
gauge symmetry breaking to Gsm- 

Acknowledgments 

We would like to thank I. Antoniadis, S. Antusch, F. Bezrukov, Wan-il Park, and M. Postma for 
useful discussions. This research was supported by the MEG and FEDER (EC) grants FPA2011-23596 
and the Generalitat Valenciana under grant PROMETEOII/2013/017. 

References 

[1] P.A.R. Ade et al. \hlCE¥2lKeck Array and Planck Collaborations], Phys. Rev. Lett. 114 (2015) 101301 [arXiv: 15 

02.00612]. 

[2] P.A.R. Ade et a/. [PtoncA: Collaboration], arXiv: 1502.02114. 

[3] M.J. Mortonson and U. Seljak, J. Cosmol. Astropart. Phys. 10 (2014) 035 [arXiv:1405.5857]. 

[4] C. Cheng, Q. G. Huang, and S. Wang, J. Cosmol. Astropart. Phys. 12 (2014) 044 [arXiv:1409.7025]; 

L. Xu, arXiv: 1409.7870. 

[5] P.A.R. Ade et a/. [BICEP2 Collaboration], Phys. Rev. Lett. 112 (2014) 241101 [arXiv: 14 03.3 985]. 

[6] R. Kallosh, A. Linde, andD. Roest,L High Energy Phys. 11 (2013) 198 [arXiv: 1311.0472]; 

R. Kallosh, A. Linde, and D. Roest, J. High Energy Phys. 08 (2014) 052 [arXiv : 1405.364 6]. 


21 



[7] J. Ellis, M. Garcia, D. Nanopoulos, and K. Olive, J. Cosmol. Astropart. Phys. 05 (2014) 037 [arXiv :1403.7518]; 

J. Ellis, M. Garcia, D. Nanopoulos, and K. Olive, J. Cosmol. Astropart. Phys. 08 (2014) 044 [arXiv: 1405.0271]. 

[8] C. Pallis, J. Cosmol. Astropart. Phys. 04 (2014) 024 [arXiv:1312.3623]; 

C. Pallis, J. Cosmol. Astropart. Phys. 08 (2014) 057 [arXiv :1403.548 6]; 

C. Pallis, y. Cosmol. Astropart. Phys. 10(2014)058 [arXiv: 1407.8522]. 

[9] C. Pallis and Q. Shaft, Phys. Rev. D 86 (2012) 023523 [arXiv: 1204.0252]; 

C. Pallis and Q. Shaft, J. Cosmol. Astropart. Phys. 03 (2015) 023 [arXiv:1412.3757]. 

[10] C. Pallis, P/i>\v. TJev. D 91 (2015) 123508 [arXiv: 1503.05887]. 

[11] M. Kawasaki, M. Yamaguchi, and T. Yanagida, Phys. Rev. Lett. 85 (2000) 3572 [hep-ph/0 00424 3]; 

P. Brax and J. Martin, Phys. Rev. D 72 (2005) 023518 [hep-th/0504168]; 

S. Antusch, K. Dutta, and PM. Kostka, Phys. Lett. B 677 (2009) 221 [arXiv: 0902.2934]; 

R. Kallosh, A. Linde, and T. Rube, Phys. Rev. D 83 (2011) 043507 [arXiv: 1011.5945]; 

T. Li, Z. Li, and D.V. Nanopoulos, J. Cosmol. Astropart. Phys. 02 (2014) 028 [arXiv:1311.6770]; 

K. Harigaya and T.T. Yanagida, Phys. Lett. B 734 (2014) 13 [arXiv: 1403.4729]; 

A. Mazumdar, T. Noumi, and M. Yamaguchi, Phys. Rev. D 90 (2014) 043519 [arXiv: 1405.3959]; 

C. Pallis and Q. Shaft, Phys. Lett. B 736 (2014) 261 [arXiv: 1405.7 645]. 

[12] D. Baumann and D. Green, J. High Energy Phys. 09 (2010) 057 [arXiv: 1004.3801]. 

[13] 1. Ben-Dayan and M.B. Einhorn, J. Cosmol. Astropart. Phys. 12 (2010) 002 [arXiv: 10 0 9.2 2 7 6]. 

[14] K. Nakayama and E. Takahashi, y. Cosmol. Astropart. Phys. 11(2010)009 [arXiv: 1008.2956]; 

K. Nakayama andE. Takahashi, J. Cosmol. Astropart. Phys. 11 (2010) 039 [arXiv:1009.3399]. 

[15] S. Antusch et al., J. High Energy Phys. 08 (2010) 100 [arXiv: 1003.3233]; 

L. Heurtier, S. Khalil, and A. Moursy, arXiv :1505.07366. 

[16] J.L. Cervantes-Cota and H. Dehnen, Phys. Rev. D 51 (1995) 395 [astro-ph/9412032]; 

M. Aral, S. Kawai, and N. Okada, Phys. Rev. D 84 (2011) 123515 [arXiv:1107.4767]; 

M.B. Einhorn and D.R.T. Jones, y. Cosmol. Astropart. Phys. 11(2012)049 [arXiv: 1207.1710]; 

J. Ellis, H.J. He, and Z.Z. Xianyu, Phys. Rev. D 91 (2015) 021302 [arXiv:1411.5537]; 

J. Ellis, T.E. Gonzalo, J. Harz, and W.C. Huang, J. Cosmol. Astropart. Phys. 03 (2015) 039 [arXiv: 1412.1460]. 

[17] M.B. Einhorn and D.R.T. Jones, J. High Energy Phys. 03 (2010) 026 [arXiv:0912.2718]; 

H.M. Lee, y. Cosmol. Astropart. Phys. OS (2010) 003 [arXiv: 1005.2735]; 

S. Ferrara et al., Phys. Rev. D 83 (2011) 025008 [arXiv:1008.2942]; 

C. Pallis andN. Toumbas, J. Cosmol. Astropart. Phys. 02 (2011) 019 [arXiv:1101.0325]. 

[18] C. Pallis and N. Toumbas, J. Cosmol. Astropart. Phys. 12 (2011) 002 [arXiv :1108.1771]; 

C. Pallis and N. Toumbas, arXiv: 1207.3730. 

[19] D.S. Salopek, J.R. Bond, and J.M. Bardeen, Phys. Rev. D 40 (1989) 1753; 

R. Fakir and W.G. Unruh, Phys. Rev. D 41 (1990) 1792. 

[20] J.L. Cervantes-Cota and H. Dehnen, Nucl. Phys. B 442 (1995) 391 [astro-ph/9505069]; 

F.L. Bezrukov and M. Shaposhnikov, Phys. Lett. B 659 (2008) 703 [arXiv: 0710.3755]; 

A.O. Barvinsky et al., J. Cosmol. Astropart. Phys. 11 (2008) 021 [arXiv :0809.2104]; 

A. De Simone, M.P. Hertzberg, and F. Wilczek, Phys. Lett. B 678 (2009) 1 [arXiv:0812.4946]. 

[21] R. Kallosh and A. Linde, J. Cosmol. Astropart. Phys. 11 (2010) 011 [arXiv:1008.3375]; 

R. Kallosh, A. Linde, and A. Westphal, Phys. Rev. D 90 (2014) 023534 [arXiv:1405.0270]. 

[22] R. Kallosh, A. Linde, and D. Roest, Phys. Rev. Lett. 112 (2014) 011303 [arXiv:1310.3950]. 

[23] S.R. Coleman and E.J. Weinberg, Phys. Rev. D 7 (1973) 1888. 

[24] K. Enqvist and M. Karciauskas, J. Cosmol. Astropart. Phys. 02 (2014) 034 [arXiv: 1312.5944]. 


22 



[25] J.L.F. Barbon and J.R. Espinosa, Phys. Rev. D 79 (2009) 081302 [arXiv:0903.0355]; 

C.P. Burgess, H.M. Lee, and M. Trott, J. High Energy Phys. 07 (2010) 007 [arXiv:1002.2730]; 

M..P.He.YX.zhsxg, J. High Energy Phys. 11 (2010) 023 [arXiv : 1002.2 995]. 

[26] F. Bezrukov, A. Magnin, M. Shaposhnikov, and S. Sibiryakov, J. High Energy Phys. 01 (2011) 016 [arXiv: 10 08 . 
5157], 

[27] A. Kehagias, A.M. Dizgah, and A. Riotto, Phys. Rev. D 89 (2014) 043527 [arXiv: 1312.1155]. 

[28] G. Lazarides and Q. Shall, Phys. Lett. B 258 (1991) 305; 

K. Kumekawa, T. Moroi, andT. Yanagida, Prog. Theor. Phys. 92 (1994) 437 [hep-ph/9405337]; 

G. Lazarides, Q. Shaft, and N.D. Vlachos, Phys. Lett. B 427 (1998) 53 [hep-ph/ 9706385]; 

G. Lazarides and N.D. Vlachos, Phys. Lett. B 459 (1999) 482 [hep-ph/9903511]; 

G. Lazarides, hep-ph/9905450. 

[29] RA.R. Ade et al. [Planck Collaboration], Astron. Astrophys. 571 (2014) A25 [arXiv :1303.5085]. 

[30] G.R. Dvali, Q. Shaft, andR.K. Schaefer, Phys. Rev. Lett. 73 (1994) 1886 [hep-ph/9406319]; 

C. Pallis and Q. Shaft, Phys. Lett. B 725 (2013) 327 [arXiv:1304.5202]; 

M. Civiletti, C. Pallis, and Q. Shaft, Phys. Lett. B 733 (2014) 276 [arXiv: 1402.6254]. 

[31] T. Dent, G. Lazarides, and R. Ruiz de Austri, Phys. Rev. D 69 (2004) 075012 [hep-ph/0312033]; 

G. Lazarides, Lect. Notes Phys. 720 (2007) 3 [hep-ph/ 0 601016]. 

[32] C. Pallis, Po5CORFU2012 (2013) 061 [arXiv: 1307 . 7 815]. 

[33] D.P George, S. Mooij, and M. Postma, J. Cosmol. Astropart. Phys. 11 (2014) 043 [arXiv: 1207.6963]; 

D. P George, S. Mooij, and M. Postma, J. Cosmol. Astropart. Phys. 02 (2014) 024 [arXiv:1310.2157]. 

[34] D.H. Lyth and A. Riotto, PAyi. Repf. 314(1999) 1 [hep-ph/9807278]; 

G. Lazarides, J. Phys. Conf. Sen 53 (2006) 528 [hep-ph/0607032]; 

A. Mazumdar and J. Rocher, Phys. Rept. 497 (2011) 85 [arXiv: 1001.0993]; 

J. Martin, C. Ringeval, and V. Vennin, Phys. Dark Univ. 5 (2014) 75 [arXiv:1303.3787]. 

[35] M.S. Turner, Phys. Rev. D 28 (1983) 1243. 

[36] D.H. Lyth, Phys. Rev. Lett. 78 (1997) 1861 [hep-ph/9606387]; 

R. Easther, W.H. Kinney, and B.A. Powell, J. Cosmol. Astropart. Phys. 08 (2006) 004 [astro-ph / 0601276]; 
D.H. Lyth, J. Cosmol. Astropart. Phys. 11 (2014) 003 [arXiv: 1403.7323]. 

[37] A. Kehagias and A. Riotto, PAyi. Rev. D 89 (2014) 101301 [arXiv: 1403.4811]. 


23 



